Disorder Induced Transitions in Layered Coulomb Gases and Superconductors 
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A 3D layered system of charges with logarithmic interaction parallel to the layers and random 
dipoles is studied via a novel variational method and an energy rationale which reproduce the known 
phase diagram for a single layer. Increasing interlayer coupling leads to successive transitions in 
which charge rods correlated in N > 1 neighboring layers are nucleated by weaker disorder. For 
layered superconductors in the limit of only magnetic interlayer coupling, the method predicts and 
locates a disorder-induced defect-unbinding transition in the flux lattice. While N = 1 charges 
dominate there, N > 1 disorder induced defect rods are predicted for multi-layer superconductors. 



Topological phase transitions induced by quenched dis- 
order are relevant for numerous physical systems. Such 
transitions are likely to shape the phase diagram of type 
II superconductors. It was proposed |lj] that the flux 
lattice (FL) remains a topologically ordered Bragg glass 
at low field, unstable to the proliferation of dislocations 
above a threshold disorder or field, providing one sce- 
nario for the controversial "second peak" line An- 
other scenario Q is based on a disorder-induced decou- 
pling transition (DT) responsible for a sharp drop in the 
FL tilt modulus. Furthermore, for the pure system, it 
was shown recently js| that in the absence of Josephson 
coupling, point "pancake" vortices, i.e vacancies and in- 
terstitials in the FL, are nucleated at a temperature Td e f, 
distinct from melting above some field. It is believed that 
this pure system topological transition merges with the 
thermal DT ||[7]] once the Josephson coupling is finite, 
being two anisotropic limits of the same transition || (at 
which superconducting order is destroyed while FL posi- 
tional correlations are maintained). Thus an interesting 
possibility is that a similar, but now disorder-induced, 
vacancy-interstitial unbinding transition can be demon- 
strated in 3D layered superconductors, relevant to many 
layered and multilayer materials (§U^]. 

In 2D recent progress was made to describe disorder in- 
duced topological transitions, in terms of Coulomb gases 
of charges with logarithmic long range interactions. It 
was shown |I(]-[l3|] that quenched random dipoles lead to 
a transition, via defect proliferation, at a finite threshold 
disorder, even at T = 0. 

In this Letter we develop a theory for a 3D defect- 
unbinding transition in presence of disorder. It is 
achieved for systems which can be mapped onto a layered 
Coulomb gas with quenched random dipoles, in which 
the interaction energy between two charges on layers n 
and n! is 2J n _ n /lnr with r the charge separation par- 
allel to the layers. One physical realization is the FL 
in layered superconductors with only magnetic 

coupling, for which we predict and locate the vacancy- 
interstitial unbinding transition. Indeed, as we argue, 
disorder induced deformations of the lattice result in ran- 
dom dipoles as seen by the defects. To study this prob- 
lem we develop an efficient variational method which al- 



lows for fugacity distributions, known |13|| to be impor- 
tant in 2D as they become broad at low T. We test 
the method on a single layer and reproduce the phase 
diagram, known from renormalization group (RG) with 
a T = disorder threshold a cr = 1/8 @. For the 
2-layer system we find that above a critical anisotropy 
rj = —J\/Jq = rj c = 1 — -j= the single layer type transi- 
tion is preempted by a transition induced by bound states 
of two pancake vortices on the two layers with a cr < 1/8. 
We develop a T = energy rationale by an approximate 
mapping to a Cayley tree problem and find that it repro- 
duces the 2-layer result. Extension to many layers with 
only nearest layer coupling shows a cascade of transitions 
in which the number of correlated charges on TV neighbor- 
ing layers increases, while the critical disorder decreases 
with rj, with N — > oo, a cr — > as i] — > 1/2. Finally we 
consider arbitrary range no for J n with the constraint 
Y) J n — 0, as appropriate for layered superconductors. 
For N > n,Q states with a cr ~ Uq/N — ► are possible 
but only at exponentially large length scales for n ^> 1. 
Thus for layered superconductors we expect that the N=l 
state dominates and find its phase diagram. Varying the 
system parameters by forming multilayers reduces no and 
allows for realization of the new N > 1 phases. 
We study the Hamiltonian: 

H = -\ Yl X! 2 ^-„'S„(r) ln(r - r')s n , (r') (1) 

(2) 



E^( r ) s «( r ) 



where s„(r) = ±1,0 define the positions r of charges 
on the n-th layer, V n (r) is a disorder potential with long 
range correlations V n (q)V n i (— q) = 47tc7Jq A„_„'/g 2 with 
A = 1 (the short distance cutoff being set to unity) . For 
simplicity we start with uncorrelated disorder from layer 
to layer A„_„/ = S nn ' with 



\V n (r)-V n (T')]2 =4aJ 2 ln|r-r'| 



(3) 



representing quenched dipoles on each layer. At T = 
the problem amounts to find minimal energy configura- 
tions of charges in a logarithmically correlated random 
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potential. For a single layer it was studied either using 
[10,[l6| a "random energy model" (REM) approximation 
[ 17 1 , or more accurately using a representation in terms of 
directed polymers on a Cayley tree (DPCT) [Q shown to 
emerge |l3| (as a continuum branching process) from the 
Coulomb gas RG of the single layer problem. Schemat- 
ically, the tree has independent random potentials (Fig. 
1) Vi on each bond with variance vf = 2a Jq . After I gen- 
erations one has 



.,2/ 



2a J$. 

sites which are mapped onto a 2D 



layer, i.e. two points separated by r ~ e have a common 
ancestor at the previous I ps In r generation. Each point r 
has a unique path on the tree (DP) with v\, ...,vi poten- 
tials and is assigned a potential V(r) — vi + ...+«;. Since 
all bonds previous to the common ancestor are identi- 
cal [V(r) — y(r')] 2 = 2Y^ i=1 Vi reproducing (||) on each 
layer. Exact solution of the DPCT |l8| yields the energy 
gained from disorder V m in — min r V(r) ps —^/8aJo\nL 
for a volume L 2 , with only 0(1) fluctuations Q, i.e 
— y8aJo per generation I = InL. 




FIG. 1. Critical disorder values with only nearest neigh- 
bor coupling Ji vs. the anisotropy r\ = — Ji/Jo. Transitions 
between different TV phases are marked with arrows. Inset: 
the Cayley tree representation (for TV = 3 neighboring layers) 
with + charges (at the tree endpoints) separated by L e along 
the layers, and separated by L from the TV = 3 — charges. 

Optimal energy configurations for M coupled layers 
are constructed considering TV neighboring layers with a 
+, — pair on each layer and no charges on the other lay- 
ers. We can take J > and J n ^o < so that equal 
charges on different layers attract. The DPCT represen- 
tation now involves, on a single tree, TV -I- polymers (each 
seeing different disorder) and TV — polymers (each seeing 
opposite disorder — Vi to their + partner). A plausible 
configuration is that the + charges bind within a scale 
L c (0 < e < 1), so do the — charges, while the + to — 
charge separations define the scale L. Its tree represen- 
tation (Fig. 1) has 2TV branches with elnL generations, 
i.e. an optimal energy of — 2N^/Sa Joe In L. On the scale 
between L e and L the 4- charges act as a single charge 
with a potential J2n=i ^«( r ) ( tne N polymers share the 
same branch) of variance TVcr hence the optimal energy is 
-2v / 8TVctJ (1 - e) lni. The total disorder energy is |§: 

E dis re -2 J V8a~[eN + (1 - e)VTV] In L . (4) 

The competing interaction energy Ei nt is the sum of 



the one for the +- pairs, [2Jo-/V+4^ =1 J n {N-n)] InL 

and for the + + / pairs, — 4Y^ =1 J„(TV — n)ehxL. 

The total energy E to t = Edi S + Ei nt being linear in e, 
its minimum is at either t = 1 or e = 0. Since e = 1 
implies that the + charges unbind, it is sufficient to con- 
sider e — with all TV > 1, i.e. a rod with TV correlated 
charges has energy (with rj n = —J n /Jo): 



Et„t = 2J N[1 - 2 



N 

E 

71=1 



Vn(l-^)-J^]\nL. (5) 



Disorder induces the TV vortex state at the critical value: 

N 



TV x - n 

°cr = "£-[1-2 2^ 



71=1 



TV' 



(6) 



(i.e. Etot — 0). Consider first only nearest neighbor cou- 
pling rji = rj\5i\. Then a cr is minimal at TV = 1 with 
a cr = 1/8 if 771 < 1 — 1/V2- For larger anisotropics suc- 
cessive TV states form at 1/(1 - 2ni) = 1 + ^N^N - 1) ~ 
TV with diverging TV as 571 -> | (Fig 1) plf. 

Consider now J„ of range no constrained by J2 n = 
as for the superconductor, e.g. r\ n — Tyie"^™ -1 ^™ for 
which a cr = (1 - e" Ar /"°)/8TV(l - e" 1 /" ). For n > 1, 
each r/n^o is small: for TV < no the lowest a cr is at TV = 1. 
However, the combined strength of TV ps no vortices be- 
ing significant a cr has a maximum and decreases back to 
zero for TV > n as a cr ps n^/8N. Hence a cr — ► as 
TV — > 00 and any small disorder seems to nucleate such 
vortices. This is because the perfect screening of the zero 
mode J2 n Jn = implies that an infinite charge rod has 
a vanishing In r interaction; hence a logarithmically cor- 
related disorder is always dominant. 

The realization of the large TV rods depends, however, 
on the type of thermodynamic limit. Adding to ([5]) the 
core energy E C N and minimizing yields a TV-vortex scale 



L w exp{E c VN/[2J Q {V8a - V8a~)}} 



(7) 



Hence as a — > such states are only achievable when 
L/N diverges exponentially. Using a cr ps tiq/8N, for 
TV > ng/8er the lowest scale L in this range is achieved 
at TV = n^/2a and leads to a lower bound L m i n ps 
cxp[75 c no/4 Jo<t] for observing large TV states with a given 
cr < |. For layered superconductors E c /J ^> 1 |^] and 
n ^> 1 and this large TV instability occurs at unattain- 
able scales, thus TV = 1 dominates. One needs no « 2 — 3, 
as in multilayers, to realize the TV > 1 states. 

To substantiate these results we develop a variational 
method for M layers which allows for fugacity distribu- 
tions, an essential feature in the one-layer problem. Dis- 
order averaging (^|) in Fourier using replicas yields: 

PHr = £pl [ s a (q,k)(GoU(q,k)st(q,k) (8) 

Jk Jq 

+(3E c J2sL(r) 
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where (G ) ab (q,k) = (4ir/q 2 )\g(k)6 ab - aJ 2 p 2 A(k)}, 
g(k) = f3J{k) — f3d^2 n J n exp(ikdn), d the interlayer 
spacing |2^] (for uncorrelated layers A(fc) = d), a,b = 
1 , . . . , m are replica indices and m — > is to be carefully 
taken. In transforming to a sine-Gordon Hamiltonian ||] 
it is crucial to keep all charge fugacitics fbl , which yields: 



(MSG = \ Xa(q,fc)(Go)^x:(q,fc) 

Jkq 



-^^y[s]exp(is • x(r)) 

r s^O 



(9) 



iFrom now on s = {s„ a }n=i : .M,a=i : .m is an integer 
vector both in layer label and replica space (i.e. of 
length m M) of entries 0,±1 and the summation is 
over all such non null vectors (also %(r) = {x«,a(r)}, 



X = Er 



iXna)- We now look for the best gaus- 



sian approximation of (^J) with propagator G } (q, k) = 
( G o)ab(^ k)+a c (k)S ab + a a (k). The bare fugacity being 
Y[s] — Gxp(—f3E c J2 na s n,a) the naive approach would 
be to restrict to charges s with a single non zero entry, 
leading to a uniform fugacity term — y E r „ a cos (x n a( r )) 
and a diagonal fc-independent replica mass term. In- 
stead we keep all composite charges s, which allow for 
variational solutions with off diagonal and /c-dependent 
replica mass terms. This corresponds respectively to 
fluctuations of fugacity and N > 1 charge rods be- 
ing generated and becoming relevant as also seen from 
RG. The variational free energy is T var — J-q + 
(Wsg — ^o)o where (...) is an average using /3Ho = 
l/ q ,fcXa(q, k)G ab (q,k)x* b ((l, k) and (3T = -±trlnG. 
The Gaussian average F[s] = y[s](expis • x( r ))o yields: 



F[s] = exp{- 



2c/ 2 



^2(G c (k)S ab -A(k))s a (k)s* b (k)} 



G c (k) = g(k) ln[A/(47r 3 (fc)a c (fc))] (10) 
A(k) = af3 2 J 2 A(k)(G c (k)/g(k) - 1) + g(k)a (k)/a c (k) 

where s a (k) = dJ2 n s na e %knd , G c (k) = G c (k) + 2^E c d, 
A the UV cutoff on q 2 . F var is minimized by a c (k)5 ab + 
fjo(fc) = kd- 2 J2 s s a{k)s* b (k)F[s\. Writing the A(k) term 
as an average over M random gaussian fugacitics Wk '■ 

exp{i|^s a (fc)| 2 A(fc)} = (exp«; fe ^s a (fc)) w (11) 



where (...)„ = JJk I -e" KI ,2Mk) d 2 w k / ^/2nA{k), 
allows to perform the exact sum on repli- 
cas yielding £ s F[s] = (Z m ) w with Z = 
E W= o, ± i } exp(-^/ fc G c (fc)| S (fc)| 2 + U>^s*(k)). 
The variational equations for m — ► become [|20f| 

. , d 2 hiZ , . , ,.dln.Z. . , , 

^ = A <ft^>- 5 ^*)= A <IW > W (12) 

For a single layer k = and Z = 1 + e u+w + e u ~ w , 
2u = —G c (0)/d, w = wq is a trinomial. (On can be 



solved for the critical line where cr c (0) — > 0. The phase 
diagram shown in Fig. 2 (full line) reproduces precisely 
recent RG results. The variational scheme, allowing for 
all replica charges s, therefore treats disorder correctly. 
For two layers kd = 0, 7r we need two fugacity distribu- 
tions wqjWk and Z is a "ninomial", i.e. Z = 1+ eight 
exponentials involving G c (0), G c (ir). Focusing on the 
low T boundary, where a c (n) ~ [er c (0)] a — > we find pc| ] 
either (i) a = 1 for 771 < 7^ c = 1 — l/y/2, representing 
decoupled layers, or (ii) a — ► 00 for 771 > 77c representing 
a ++ bound states on the two layers. The T = energy 
rationale is therefore reproduced. The phase diagram for 
two layers with r\ c < 77 < 1/2 is shown in Fig. 2 fl5|| 




T/J 

FIG. 2. Phase diagram for the onset of the N = 1, 2 insta- 
bilities for anisotropy 77 = 0.35. At low T two distinct transi- 
tions are possible, the first being to the rod N = 2 phase. At 
high T the independent layer N — 1 transition dominates 

For any number of layers one obtains a simple N rod 
solution by restricting the sum over s in @ to a sub- 
class of charges of the form s na = s a J2j=i,N <*n,n'+j-i- 
The variational solution, of the form u c (k) = <r c 4>N{k), 
reduces to an effective one layer problem, in term of 
the structure factor of the rod s a (k)sl(k) — <j>N(k) = 
sin 2 (Affcd/2)/sin 2 (£;d/2). The N rod becomes critical at: 



(?cr = ( / <Mfcy(fc))7(8J / &v(*)A(fc)) (13) 

a formula which can equivalently be obtained within the 
Cayley tree rationale. Indeed for any correlations A„, 
the energy of the e = configurations is still given by (^) 
replacing a^a(l + 2 £^ =1 A„(l - n/N))). © repro- 
duces both single (using f/>(0) = iV) and two layer results. 
Finally, for many layers and weak interlayer coupling (e.g. 
771 < 1 — l/v2 in Fig.l) the N = 1 transition dominates 
and occurs at a cr = 1/8 ((|T^) using ^i(fc) = 1). 

As a direct application we consider a flux lattice 
in a layered superconductor with no Josephson cou- 
pling and a magnetic field B perpendicular to the lay- 
ers. The FL is composed of pancake vortices displaced 
from the p-th line position R p at the n-th layer into 
R p + Up. The defects s n (r) couple to the lattice via 
H-vac = J2r P nn> s„> (r)G v (R p + - r, n' - n) where, in 
Fourier § G£(q,fc) = (0 2 d 2 /4^A 2 bt2 2 )/[l + /(q, k)} where 
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/(q,fc) = (d/4\l b q)smhqd/[sinh 2 (qd/2) +sin 2 (kd/2)}-, 
cba = Ba 2 is the flux quantum, a the FL spacing, A a fc 
the penetration length along the layers. To 0-th order in 
u™ the defects feel a periodic potential fixing their posi- 
tion in a unit cell, hence s(q, k) involve only \q\ < 1/a. 

In the limit q — ► the longitudinal modes, to 
which defects couple, have for (tilt) elastic energy 

11 W ei = 53^ / fc ,-D(*0|UL(cbfc)| a With D(k) = 

|Eq/o[ G *(Q> fe ) - Gv(Q,0)] + G v (k) where Q arc 
reciprocal wavevectors of the lattice and G v (k) = 
lim q ^ Q G v (q,k)q 2 = ^d 2 fc 2 /[47r(l + A 2 b fc 2 )] and k z = 
(2/d) sin(fcd/2). The sum on Q is due to the high mo- 
mentum components of the magnetic field and is respon- 
sible for the non-perfect screening of the defect interac- 
tion and to a finite Td e f ■ Minimizing H va c + 7~L e i yields 
Ut, ac (q, k) = «qs(q, k)G v {k)a 2 / D(k)q 2 and (@) with: 



g(k) = (3G v (k)[l - G v (k)/D(k)}/4Tr 



(14) 



g(k) = (3dr' ln(l + a 2 fc 2 /47r) with r' 



and | T t 



def 
2, 



Thus the long range interaction is ~ lnr and its coef- 
ficient determines Td e f = 2 Jo (via J k g(k) = 2). Since 
fj. G v (k) ~ <l>od/\ 2 b , the scale of the melting transition 
|L4|, the defect transition occurs before melting and can 
thus be consistently described only if D(k) — G v (k) -C 
D(k). This is possible if either d <C a <C A a {, where 

2 da 2 /(128^A 4 ab ) 
= T r ln(a/d) or for d > a where g(fc) = 
l3{d 4 /a)T'k 2 z e- 2 * d / a leading to T de/ = 4(d/a)r' e - 27rd / a . 
Remarkably -D(fc) — G v (k) <C -D(A:) also yields that the 
long range response u„ ac (r) ~ a 2 r/r 2 to a vacancy at 
r = is confined to the same layer. 

Point disorder deforms the flux lattice, producing 
quenched dipoles coupling to our defects. Expansion of 
the disorder energy, valid below the Larkin length |Q, 
and minimization together with 7i vac + 7~L e i yields read- 
ily (|^). A more general argument, valid at all scales, 
treats u vac as a small perturbation around the Bragg 
glass configuration. Systematic expansion of the free en- 
ergy F = F B g + gajjr J qk iqs(q, k)G v (q, k)(u(q, k)) s=0 + 
0(s 2 ) in defect density in a given disorder config- 
uration shows that a defect feels a logarithmically 
correlated random potential Vi(r) as in (0, Isl) with 
aJ^A(k) = G v (k) 2 lim q ^ a C B G(q,k)/4Trd 2 a^ where 
Cbg{ v iI) = (uq (0)) (uf (r)) is the correlation in the un- 
perturbed Bragg glass Cbg(0, k) ~ l/(c 2 4 (fc 4 + i?^ 1 /? 3 ), 
R c a Larkin length along c iQ . It yields a k- independent 
A(fc) for k>l/R c while A(fc) ~ k for k < l/R c . 

Applications to FL depends on the interlayer form 
of (|l4| ) of range no ~ a/d for large a/d. Remarkably 
g(k = 0) = 0, i.e perfect screening holds as in 2D ||. 
Hence ^ n J n = and as uq is reduced Jo, J\ dominate 
the sum, i.e. t]\ — > \ when d 3> a. One finds that rji 
crosses the critical value 1 — l/v2 when d/a » 1, de- 
pending weakly on a/\ a b- We thus propose that FL in 
multilayer superconductors, where d > a can be achieved, 
can show a rich phase diagram with N > 1 phases. 



In layered superconductors a/d w 10 — 100 Q and the 
N = 1 transition at a cr — 1/8 dominates for realistic 
sizes. The disorder-induced decoupling transition, ne- 
glecting defects, predicted at cr^ec = 2 is thus above 
the defect transition (with B ~ a) in the B — T plane 
(similarly thermal decoupling occurs at Td ec = 8Td e f for 
d<a«A). A natural scenario is again of a single tran- 
sition at <j c varying from 2 to 1/8 as the bare Josephson 
coupling is reduced, e.g. by increasing d in multilayers. 

In conclusion, we developed a variational method and 
a Cayley tree rationale to study layered Coulomb gas. 
The results are relevant to flux lattices where we find 
the phase boundaries and propose new N > 1 phases for 
d> a. The present methods may be useful for other 2D 
disordered systems, such as quantum Hall. 
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